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Surface velocity network with anti-solar differential rotation on the ac- 
tive K-giant a Geminorum 
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We demonstrate the power of the local correlation tracking technique on stellar data for the first time. We recover the spot 
migration pattern of the long-period RS CVn-type binary a Gem from a set of six Doppler images from 3.6 consecutive 
rotation cycles. The resulting surface flow map suggests a weak anti-solar differential rotation with a ~ —0.0022 ± 
0.0016, and a coherent poleward spot migration with an average velocity of 220 ± 10ms" 1 . This result agrees with our 
recent findings from another study and could also be confirmed theoretically. 

© 2007 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



1 Time-series Doppler images of a Gem 



a Geminorum (75 Gem, HR2973, HD 62044) is a long- 
period (P ro t ~ 19.6 days) RS CVn-type system with a Kl- 
giant primary and an unseen companion. In our initial paper 
(Kovari et al. 2001) six time-series Doppler images were 
recovered from the Ca I-6439A and Fe I-6430A lines. The 
maps cover 3.6 consecutive stellar rotations in 1996/97. Im- 
ages #1, #3, and #5, as well as #2, #4, and #6 represent 
contiguous stellar rotations and therefore independent maps 
(see our extensive study for details). Since there is a rea- 
sonable good agreement between the respective Ca and Fe 
maps, for further use we combine them to decrease statisti- 
cally the uncertainties of the individual Doppler reconstruc- 
tions. 

2 Local correlation tracking 

The local correlation tracking (hereafter LCT) technique, 
originally developed for tracking different solar surface fea- 
tures, cf., e.g., Sobotka et al. (1999, 2000), Ambroz (2001), 
Svanda et al. (2006), is based on the principle of the best 
match of two image frames that record the tracked surface 
structures at two different instants. A small correlation win- 
dow is chosen in the first frame and compared with a some- 
what displaced window in the second frame. When the best 
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matching displacement is found, a vector of displacement 
is recorded, and a new correlation window is chosen by the 
side of the former one. This pattern is followed until the 
whole image frame is covered. The resulting vector field 
is interpreted as the surface flow map over the time lag of 
the two initial images. For a more detailed description of the 
method as well as numerical simulations with artificial data, 
we refer to a forthcoming paper (Svanda et al. 2008). 

In our pilot study we apply the LCT technique to the 
time-series of the six average (Ca+Fe) Doppler images of 
a Gem (Kovari et al. 2001). From the six maps we set 
four correlation pairs: #l-#3, #2-#4, #3-#5 and #4-#6, 
i.e., we choose the independent consecutive but contiguous 
maps to minimize the time shift and, thus, the masking ef- 
fect of individual spot evolution on the flow pattern. Then, 
for each pair a velocity field map is computed and finally 
the four maps are averaged. 

The resulting average flow map is shown in Fig.Q] In the 
figure the large-scale flow network is related to the spatially 
resolved surface structures, which is expected from solar 
analogy. The amplitude of the velocity vectors of typically 
several hundreds of ms _1 also seem reasonable. Flows of 
the order of 100 ms -1 are observed in the photosphere of 
the Sun. Fig. Q] shows convergent flows towards spotted re- 
gions, which is also consistent with the solar case. 

In Fig. |2] we plot the zonal and meridional flow compo- 
nents extracted from the LCT map. Despite its large error 
bars, the latitudinal distribution of the zonal component is 
in favor of a weak anti-solar differential rotation (DR after- 
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Fig. 1 Surface changes in the photosphere of a Gem from LCT. As a background we use the average of the six combined 
(Ca+Fe) time-series Doppler images. Two distinct flow structures can be observed: a complex network of convergent flows 
around spots and a general poleward trend. 



wards). The latitude (J3) dependent rotation law is assumed 
in the usual quadratic form of 

= n cq - Aft sin 2 /? , (1) 

where f2 oq is the equatorial angular velocity, Aft — J7 cq — 
fipoie is the angular velocity difference between the equa- 
tor and the pole. The surface shear parameter is defined as 
a = Aft/n cq . The best fit gives n cq = 18.33 ± 0.017day 
andAfi = -0.04± 0.037day which yields a = -0.0022± 
0.0016, i.e., an almost rigid body rotation. The meridional 
flow component indicates mostly poleward spot migration. 
The formal surface-averaged poleward meridional compo- 
nent is 2.5 ± 0.1° over a rotation cycle of w20 days, which 
converts to an average velocity of « 220 ± 10ms _1 . In the 
RS CVn-type a Gem such a meridional flow could be main- 
tained by large-scale thermal inhomogeneities (from large 
cool spots) and/or by tidal effects from the close binary na- 
ture (cf. Kitchatinov & Rudiger 2004). 



we make a linear normalization, since the time baseline of 
the ccf maps are different. After averaging the normalized 
maps, we fit the correlation peak for each latitudinal stripe 
with a Gaussian profile. These Gaussian peaks per latitude 
are fitted with the usual quadratic form (Eq. 1). The result 
shown in the upper panel of Fig. 3 suggests an anti-solar 
type DR law with fi oq = 18.29 ± 0.057day and Afl = 
-0.38 ± 0.087day and a surface shear of a « -0.021 ± 
0.005. 

Meridional motion of surface features can be derived 
similarly, by cross-correlating the corresponding longitude 
stripes along the meridian circles (but restricting only to the 
more reliable visible hemisphere). The resulting average lat- 
itudinal ccf map is plotted in the lower panel of Fig. 3. The 
best correlating latitudinal shifts suggest a common pole- 
ward migration of spots of « 4.1 ± 0.3° per rotation cycle, 
which can be interpreted as a poleward meridional flow with 
an average velocity of w 350m s _1 . 



3 Cross-correlation results 

In this section we compare the LCT results with our recent 
findings (Kovari et al. 2007b) for the same data by apply- 
ing the method of Average cross-correlation of contiguous 
Doppler images' (hereafter ACCORD). The method was de- 
scribed and applied to detect surface differential rotation 
first for LQ Hya in Kovari et al. (2004), and further on has 
become a powerful tool in analysing time-series Doppler 
maps (e.g., Kovari, Weber & Strassmeier 2005, Kovari et 
al. 2007a). 

We use again the consecutive but contiguous correlation 
image pairs (i.e., #l-#3, #2-#4, #3-#5 and #4-#6) and 
compute four cross-correlation function (ccf) maps. Then 



4 Summary and conclusions 

We have demonstrated that LCT technique can be used to 
map flow fields on stellar surfaces. Applying LCT for time- 
series Doppler images of a Gem we reconstruct large-scale 
surface flows in the order of several hundreds of ms" 1 and 
observe complex network of convergent flows around spots. 
Zonal flow components suggest rigid-body rotation or anti- 
solar DR with a very weak surface shear of a = —0.0022 ± 
0.0016, while meridional flow components show a coherent 
poleward spot migration at an average velocity of 220 ± 
10 ms -1 . For a more detailed description of the method as 
well as for numerical tests on artificial data, we refer to a 
forthcoming paper (Svanda et al. 2008). 
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Fig. 2 Zonal (top) and meridional (bottom) flow compo- 
nents from the LCT method. The zonal distribution may 
be indicative of anti-solar differential rotation with a shear 
of a — —0.0022 ± 0.0016. The meridional component 
shows a coherent poleward flow with an average velocity 
of 2.5 ± 0.1° over a rotation cycle. 

Our findings agree well with the results from a re- 
cent application of the ACCORD method on the same data 
(Kovari et al. 2007b). ACCORD resulted also in anti-solar 
DR law, but with a stronger shear of a = -0.021 ± 0.005. 
From the latitudinal cross-correlation a coherent poleward 
flow was derived with a rate of w 350 ms -1 , which is 
in qualitative and also in quantitative agreement with the 
meridional flow component from LCT. The deviation could 
be originated from the different limitations of the two dif- 
ferent approaches. 

Sufficiently fast meridional flow can result in anti-solar 
type DR and the necessary angular momentum transport 
(i.e., the necessary rate of the meridional flow) can be esti- 
mated from Eq. 15 in Kitchatinov & Riidiger (2004). Taking 
a turnover time of r = 5.5 x 10 6 s with a mixing length I 
of 7 x 10 8 m (cf. Gunn, Mitrou & Doyle 1998, Paterno et 
al. 2002) would yield « 300 ms -1 , which agrees with our 
results derived from two different methods. If the poleward 
surface migration on a Gem is interpreted due to an under- 
lying meridional circulation, our result could be regarded as 
a confirmation of the theoretical estimates in Kitchatinov & 
Riidiger (2004). 
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Fig. 3 Latitudinal (top) and longitudinal (bottom) cross- 
correlation maps of a Gem from ACCORD using the av- 
erage (Ca+Fe) Doppler images from Kovari et al. (2001, 
2007b), suggesting anti-solar DR with a coherent poleward 
migration of spots. 
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